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Abstract 



We study the late time evolution of a class of exact anisotropic cosmological so- 
lutions of Einstein's equations, namely spatially homogeneous cosmologies of Bianchi 
type VIIo with a perfect fluid source. We show that, in contrast to models of Bianchi 
type Vllh which are asymptotically self-similar at late times, Bianchi VIIo models un- 
dergo a complicated type of self-similarity breaking. This symmetry breaking affects 
the late time isotropization that occurs in these models in a significant way: if the 
equation of state parameter 7 satisfies 7 < | the models isotropize as regards the 
shear but not as regards the Weyl curvature. Indeed these models exhibit a new dy- 
namical feature that we refer to as Weyl curvature dominance: the Weyl curvature 
dominates the dynamics at late times. By viewing the evolution from a dynamical 
systems perspective we show that, despite the special nature of the class of models 
under consideration, this behaviour has implications for more general models. 

1 Introduction 

In cosmology there are two asymptotic regimes, namely 

i) the approach to the initial singularity, characterized by I — > 0, and 

ii) the late time evolution, characterized by / — > 00, 

where I is the overall length scale. In other words, the asymptotic regimes correspond to 
the extreme values of the length scale variable. Alternatively, one can characterize the 
asymptotic regimes in terms of the Hubble variable H, which is related to / by 



where the overdot denotes differentiation along the fundamental congruence. For an ever- 
expanding universe, H > throughout the evolution, and the singular regime is characterized 
by H — > +00 and the late-time regime by H — > 0. 

The two asymptotic regimes are of course totally different from a physical point of view, 
since at the singularity physical quantities (for example, the matter density, the shear of the 
Hubble flow, the Weyl curvature) typically diverge, while at late times physical quantities 
typically tend to zero. Despite this physical distinction, the asymptotic regimes of the 
simplest ever-expanding models, the isotropic Friedmann-Lemaitre models and the anistropic 
but spatially homogeneous models of Bianchi type I, have an important property in common: 
they are asymptotically self- similar [], that is, they are approximated by self-similar models 
in the asymptotic regime, as was first pointed out by Eardley (1974, page 304). 

The above behaviour becomes less surprising in view of the following heuristic considera- 
tion: in both asymptotic regimes, one has a self-gravitating system evolving in size through 
many orders of magnitude (evolution into the past for the singular regimes and into the 

1 The flat FL model is in fact (exactly) self-similar (Eardley 1974, page 304). 
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future for the late time regime), making it plausible that the system might become scale- 
invariant i.e. self-similar, asymptotically. It is thus tempting to speculate that a "Principle 
of Asymptotic Self- Similarity" holds quite generally for Bianchi universes, namely that any 
Bianchi universe is approximated by a self-similar Bianchi universe in the asymptotic regimes. 
However, the well-known Mixmaster universes (Bianchi types VIII & IX) provide a decisive 
counter-example to the validity of such a principle for the singular regime, since they oscil- 
late indefinitely as the initial singularity is approached into the past, and thus do not have a 
well-defined asymptotic state. Nevertheless the Mixmaster universes reinforce the idea that 
self-similar models are important as regards the overall evolution of Bianchi universes, since 
they evolve through an infinite sequence of Kasner (i.e. self-similar) states as they approach 
the initial singularity into the past. Indeed, research to date strongly suggests that a non- 
tilted perfect fluid Bianchi universe either approaches a unique self-similar asymptotic state, 
or it passes through an infinite sequence of self-similar states as one follows the evolution 
into the past towards the initial singularity. 

The question remains as to the validity of the Principle of Asymptotic Self-Similarity in 
the late-time regime. Our main goal in this paper is to show that the non-tilted Bianchi VIIq 
universes provide a counter-example to asymptotic self- similarity in the late time regime. 

One can gain a deeper understanding of asymptotic self-similarity in Bianchi universes 
and of the reasons for its violation by viewing the Einstein field equations from a dynamical 
systems perspective. In studying physical phenomena, it is usually desirable to use dimen- 
sionless variables. In cosmology the density parameter, a dimensionless quantity defined by 



where /i is the matter density, G is the gravitational constant and a subscript zero means 
evaluation at the present time, is one of the fundamental observational parameters. In 
addition the dimensionless quantity (o~/H)q, where a is the shear scalar of the Hubble flow, 
gives a measure of the anisotropy in the Hubble flow. This quantity is constrained by 
the observations of the anisotropy in the temperature of the cosmic microwave background 
radiation (CMBR). 

In view of these considerations, it is natural to formulate the Einstein field equations 
using dimensionless variables that are defined by normalizing in an appropriate way with 
the Hubble scalar H, as in Wainwright and Hsu 1989. The appropriate time variable is a 
dimensionless quantity r that is related to the length scale I by 



where Iq is a constant. We shall refer to r as the dimensionless time variable. In terms of r, 
the singular regime is defined by r — > — oo and the late time regime for an ever-expanding 
model by r — > +00. The dimensionless time r is related to clock time t by 









I = l e T , 




dt _ 1 
dr = H' 




2 From now on we will use geometrized units, i. e. 8ttG = 1 and c = 1 
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as follows from equations ( |1.1|) and ( |1.3| ). We refer to Wainwright and Ellis 1997 (Chapters 
5 and 6) for full details []. 

For Bianchi universes the Einstein field equations with a perfect fluid source and an equa- 
tion of state p = (7 — reduce to an autonomous system of ordinary differential equations, 
thereby defining a dynamical system. It turns out that if one uses the above dimensionless 
variables, the equilibrium points (i.e. fixed points) of this dynamical system correspond to 
self-similar Bianchi universes, provided that 7 > 0. Within this framework, the evolution of 
a Bianchi universe is described by an orbit of the dynamical system, and asymptotic self- 
similarity into the past /future means that the orbit approaches an equilibrium point into 
the past/future. Conversely, asymptotic self- similarity is violated whenever the orbit of a 
cosmological model does not approach a single equilibrium point into the past/future. 

In the Mixmaster model the asymptotic self-similarity breaking is due to the fact that 
into the past, the orbits approach a bounded two-dimensional attractor containing the set of 
Kasner equilibrium points (the so-called Mixmaster attractor, see WE, page 146), instead of 
a single equilibrium point. On the other hand, we will show that the breaking of asymptotic 
self-similarity in the late time regime in Bianchi VIIo models is due to the fact that state 
space is unbounded, and the orbits "escape to infinity" at late times. The difference between 
the two types of symmetry breaking is reflected in the behaviour of dimensionless scalars 
formed by normalizing with the Hubble scalar H. In a self-similar model such scalars are 
constant. It follows that in an asymptotically self-similar regime the limits of such scalars 
exist. In a Mixmaster universe the limits of some scalars of this type do not exist in the 
singular asymptotic regime (see WE, figures 11.1-11.3), but all remain bounded. On the 
other hand, in a Bianchi VIIo late time regime we will show that, depending on the equation 
of state of the perfect fluid, a dimensionless scalar constructed from the Weyl curvature 
tensor oscillates and becomes unbounded. 

In section 2 we introduce the three primary dimensionless scalars and summarize their 
behaviour in the late time asymptotic regime. This section contains no technical details, 
and is intended to give an overview of the main results. In Section 3 we introduce the basic 
dimensionless variables, which are based on the orthonormal frame formalism, referring to 
WE (chapters 5 and 6) for full details. A change of variable then leads to a new form of 
the evolution equations that is adapted to the oscillatory nature of the Bianchi VIIo models. 
The detailed asymptotic form of the solutions of these evolution equations in the late time 
regime is given (see theorems 3.1 - 3.3). Section 4 contains a discussion of the cosmological 
implications of our results. 

There are three Appendices. Appendix A contains a proof of the fact that Bianchi VIIo 
universes are not asymptotically self-similar at late times, while Appendix B contains the 
lengthy proof of Theorem 3.1, which gives the asymptotic form of the solutions when the 
equation of state parameter satisfies | < 7 < |. In Appendix C we use the asymptotic 
expansions in Theorems 3.1 - 3.3 to determine the asymptotic form of the line-element, so as 
to be able to compare our results with those of other investigations. We note, however, that 
knowledge of the orthonormal frame components is sufficient for analyzing any aspect of the 
model, for example the degree of isotropization and the anisotropy of the cosmic microwave 
background radiation. The asymptotic form of the line-element is thus not central to our 

3 This work will henceforth be referred to as WE. 
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conclusions, and hence is placed in an appendix. 



2 Main results 

Our main results concern the dynamics in the late time regime of non-tilted perfect fluid 
Bianchi universes of group type VII . The Bianchi VII universes that are locally rotationally 
symmetric Q either admit a group of Bianchi type I or are flat FL universes, and hence are 
asymptotically self-similar. We thus exclude this special class in the sequel. The first result 
concerns the breaking of asymptotic self- similarity. 

Theorem 2.1: 

If | < 7 < 2, then any non-tilted perfect fluid Bianchi universe of group type VIIo that is 
not LRS is not asymptotically self-similar as r —>■ +oo. 

The proof of this Theorem, which requires some technical results from the theory of 
dynamical systems, is given in Appendix A. 

As mentioned in the Introduction, the Hubble scalar H is used to define dimensionless 
variables in cosmology. We recall that H is defined by H = \u a a , where u is the 4-velocity 
of the Hubble flow (the cosmo logical fluid), and determines the lenth scale I according to 
( |1.1| ). If a Bianchi universe is asymptotically self-similar at late times, then all dimensionless 
scalars formed by normalizing with the Hubble scalar H have limits as r — > +oo. As a 
result of Theorem 2.1, it is possible that the limits of some dimensionless scalars will not 
exist as r — > +oo. The next results give the limiting behaviour of three physically important 
dimensionless scalars. 

The first scalar is the density parameter Q, defined by Q 



n=3^. (2.D 

where /i is the matter density of the perfect fluid. The density parameter gives a measure of 
the dynamical significance of the matter content of the model. If Q <C 1 then the dynamics 
will be close to a vacuum model, characterized by Q = 0. 
The second scalar is the shear parameter E defined by 

S 2 = (2.2) 

where a 2 = \cr a i J iJ ab , and a a b is the rate-of-shear tensor of the Hubble flow. The shear 
parameter gives a dimensionless measure of the anisotropy in the Hubble flow by comparing 
the shear scalar a to the overall rate of expansion as described by H. The anisotropy in the 
temperature of the CMBR enables one to estimate the value of E at the present epoch (see 
for example Maartens et. al. 1996). 



See for example WE, pages 22-3; we use the standard abbreviation LRS. 
5 It is customary to use a subscript zero on scalar such as f2 and H, as in equation Jl.2| ), to indicate their 
present day value. 
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The third scalar gives a dimensionless measure of the Weyl curvature tensor, and is 
defined by 



W 

W=Jp, (2-3) 

where 

W 2 = l(E ab E ab + H ab H ab ). (2.4) 
Here E ab and H a b are the electric and magnetic parts of the Weyl tensor, defined by 

E a b = CarbsU U , H a b = C ar b s U U S , 



where u is the 4- velocity of the Hubble flow and * denotes the dual operation. We shall refer 
to W as the Weyl parameter. W can be regarded as describing the intrinsic anisotropy in 
the gravitational field. Cosmological observations can in principle give an upper bound on 
W, although obtaining a strong bound is beyond the reach of present day observations. We 
refer specifically to observations of the CMBRQ (see Maartens et al 1995b), and observations 
of distant galaxies (see Kristian & Sachs 1966, page 398). 
The generalized Friedmann equation 

3H 2 = a 2 - \ 3 R + fi, (2.5) 

where 3 R is the scalar curvature of the spacelike hypersurfaces of homogeneity, places bound 
on n and £ 2 . Dividing by 3H 2 yields 

i = y? + k + n, (2.6) 

where 

K ■ '* 



6H 2 

is a dimensionless quantity. Since K > for Bianchi type VIIo geometries Q, equation (|2] 
implies that 

S 2 < 1 and Q < 1. (2.7) 

On the other hand there is no general upper bound for the Weyl parameter W. Indeed we 
shall show that W is unbounded in the state space of the Bianchi VIIo cosmologies. 

We now describe in turn the asymptotic behaviour of the three primary expansion- 
normalized scalars Q, £ and W in Bianchi VIIo universes at late times, i.e. as r — >• +oo. 
These results will follow immediately from Theorems 3.1 - 3.3 in Section 3. 



6 The analysis of Maartens et al 1995b does lead to a strong upper bound on the Weyl parameter W (see 
their equation (41)). However, they make an assumption concerning the time derivatives of the temperature 
harmonics (assumption C2') which is not observationally testable at present. 

7 This result can be established using equations (1.95), (1.103) and (1.105) in WE. 
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Theorem 2.2: 

For any non-tilted perfect fluid Bianchi VII universe that is not LRS, the density parameter 
Q satisfies 



. u 2 — 4 

lim f2 



if 1< 7 < 



3 " ' — 3 

'1(2-7), if |<7<2. 



Theorem 2.3: 

For any non-tilted perfect fluid Bianchi VII universe that is not LRS, the shear parameter 
E satisfies 

lim £ 2 = 0, if | < 7 < I 

and 

lim sup £ 2 = 1(37-4), liminf S 2 = 1(3 7 - 4) 2 , if \ < 7 < 2. 
Theorem 2.4: 

For any non-tilted perfect fluid Bianchi VII universe which is not LRS, the Weyl curvature 
parameter W satisfies 



lim W 

t^+oo 



, if |< 7 <1, 
L^O, if 7 = 1, 
+00 , if 1 < 7 < 2. 



where L is a constant that depends on the initial conditions, and can have any positive real 
value. 

The implications of Theorems 2.3 and 2.4 are as follows: 

i) there is a shear bifurcation at 7 = |; specifically, if 7 > | the models no longer 

isotropize as regards the shear of the Hubble flow, and in fact lim S 2 does not exist 

t^+oo 

since S 2 is oscillatory as r — > +00. 

ii) there is a PFey/ curvature bifurcation at 7 = 1; specifically, if 7 = 1, the models no 
longer isotropize as regards the Weyl curvature, and if 1 < 7 < 2 the Weyl curvature 
moreover dominates the rate of expansion at late times. 



Comment: By Theorems 2.2 - 2.4, the limits as r — > +00 of Q, S and W exist if | < 
7 < 1. One might thus expect that the models are asymptotically self-similar in this case, 
contradicting Theorem 2.1. We will show in Section 3, however, that if | < 7 < 1, a 
dimensionless scalar representing a time derivative of the Weyl tensor is unbounded as 

T — > +OO. 
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3 Evolution equations and asymptotic behaviour 



In order to write the evolution equations for the non-tilted perfect fluid Bianchi universes 
of group type VII in a suitable form, we use the ort honor mal- frame formalism, introduced 
in cosmology by Ellis & MacCallum 1969, in which the commutation functions 7^ of the 
orthonormal frame {e a }, defined by 

t e a, e 6 ] = 7^e c , 

act as the basic variables of the gravitational field. We choose the frame to be invariant 
under the isometry group G3, and aligned with the fluid velocity u in the sense that e =u. 
It follows that the 7^ depend only on t, the clock time along the fluid congruence. For 
Bianchi universes of group type VIIo one can specialize the frame so that the only non-zero 
commutation functions are the Hubble scalar H = |6, where 9 is the rate of expansion 
scalar, and the diagonal components of the shear tensor <j a p and a matrix n a p, 

a al3 = diag((7n, 0-22, 0-33), n afS = diag(0, n 2 , n 3 ), (3.1) 

with n 2 > and n 3 > 0. (see Ellis & MacCallum 1969). Since a a p is trace-free, it has only 
two independent components, which it is convenient to label using 

0"+ = |(o"22 + C33), C- = iT^ " 22 ~~ °" 33 )- ( 3 - 2 ) 

In analogy we label the components n 2 ,n 3 by 

n + = ±(n 2 + n 3 ), n_ = ^(n 2 - n 3 ) (3.3) 
The physical state of the models is thus described by the variables 

(H,a + ,a_,n + ,n_). (3.4) 
We now introduce expansion-normalized variables and a dimensionless time r according 

to 

E ± = f, = (3.5) 

and replace t by the dimensionless time variable r according to equation (1.4) (see WE, 
Chapter 5, for background and motivation). The evolution equation for H can now be 
written as 

H' = -(l + q)H, (3.6) 

where ' denotes differentiation with respect to r, and q is the deceleration parameter, given 
by 

g = 2£ 2 + i(3 7 -2)a (3.7) 
(see WE, pages 113-4). In addition S 2 and Q are given by 

S 2 = S 2 F + S 2 _, (3.8) 
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= 1 - l£ - £i - JVi. (3.9) 



The evolution equations for E-±. and A^ are 

E' + = -(2 - g)E + - 2A^ 
E'_ = -(2-g)E_ -2N-N+ 

iV; = (g + 2£ + )iV + + 6E_AL (3.10) 
iVi = (g + 2E+)AL + 2£_iV + . 

These equations follow from equations (6.9), (6.10) and (6.35) in WE on setting Aq = 0. For 
future reference we also note the evolution equation for Q: 

Q' = [2q - (37 - 2)]fl, (3.11) 

(see WE, page 115). 

The dimensionless state vector is (E+, E_, A^ + , iV_) and the physical region of state space 
is defined by the condition Q > 0, i.e. 

T,\ + T?_ + Ni < 1, (3.12) 

and the restrictions 

Ni-3N 2 >0, N + >0, (3.13) 



which are a consequence of ( |3.3|) , (|3.5|) and the restrictions n 2 > 0, n 3 > 0. The variables 
E_|_, E_ and AL are thus bounded, but AT+ can take on any real value. Indeed, we show in 
Appendix A that if Q > and | < 7 < 2, then for any initial conditions, 

lim N + = +00. (3.14) 

t^+oo 

The variable N + does not appear in the equation for E' + , but does play a significant role 
in the equations for E'_ and N'_, which are of the form 



E'_ = —2N^N + + {bounded terms} 
N'_ = 2E_A^ + + {bounded terms} . 

The leading terms generate oscillations in E_ and AL, and suggest that we introduce 
polar coordinates in the E_AL-space. We also replace the variable AL by 1/AL so as to 
have a variable that remains bounded as r — > +00. We thus let 

E_=#cos^, AL = .Rsin^, (3.15) 

where R > 0, and 




(3.16) 
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where 



A straightforward calculation using (3.10), ( |3.15|) and ( |3.16|) yields the evolution equa- 
tions for the new variables (E + , R, M, ip) in the following form: 

E' + = -(2-Q)£ + - J R 2 + (l + £ + ) J R 2 cos2?/; (3.17) 

R' = [Q + £+ - 1 + (R 2 - 1 - £+) cos 2if>] R (3.18) 

M' = -[Q + 2£ + + i? 2 (cos 2if) + 3M sin 2^)] M (3.19) 

^ = j-[2 + (l + S + )M sin 2^] , (3.20) 

Q = 2£ 2 ^ + J R 2 + i(3 7 -2)ft, (3.21) 

= 1 - - R 2 (3.22) 

In writing these equations we have expressed the trigonometric dependence in terms of 
cos 2if> and sin 2ip. In particular, we have decomposed the deceleration parameter q, as given 
by equation (|3.7|), into a non-oscillatory and an oscillatory part, 

q = g + J R 2 cos2^. (3.23) 

The variables £+, R, M and if) are required to satisfy the following restrictions. First, the 
requirement Q > is equivalent to 

S^ + i? 2 <l. (3.24) 
Second, the restrictions ( |3.13|) in conjunction with the definitions ( 3.15|) and ( 3.16 ) lead to 



M > 0, R > 0, 3R 2 M 2 sin 2 ip < 1. (3.25) 

We note that if R = equation (3.20) for if}' becomes irrelevant, and the remaining evolution 
equations (3.17) and (3.19) describe the LRS Bianchi VII models (which are equivalent to 
LRS Bianchi I models). 

We now state the main results concerning the asymptotic behaviour of the solutions of 
the DE (3.17) - (3.20) with R > 0, Q > and | < 7 < 2. For convenience, we introduce 
a new parameter f3 given by 

0= 1(4-37), (3-26) 

and note the relations 

§(2- 7 ) = 1+0, i(37-2) = l-0, 3( 7 -l) = l-20. (3.27) 
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Theorem 3.1: 

!<7<f 



If | < 7 < f , any solution of the DE (3.17) - (3.20) with R > and Q > satisfies 



lim M = 0, lim R = 0, lim £_ 



T— >+0O 



(3.28) 



and 



M = C M e~^ 



R = Cne-P 



c -2/3t 



1-/3 



l + 0(e" fer ) 

i + o( e - fer ) 

l + 0(e- 6r ) 



(3.29) 



as r — ► +oo, where /3 is given by ( |3.26| ), Cr and Cm are constants that depend on the initial 
conditions, and h is a positive constant. 



Theorem 3.2: 

4 
3' 



If 7 = |, any solution of the DE (3.17) - (3.20) with R > and SI > satisfies 



M = C M re- T 
R- 1 



r 2 



2 ' 



i + o(^) 

l + 0(^i) 

i + o(^) 



^81), and 



(3.30) 



where Cm is a constant that depends on the initial conditions. 



Theorem 3.3: 

If | < 7 < 2, any solution of the DE (3.17) - (3.20) with R > and Q > satisfies 



lim M = 



lim £_ 



T— >+0O 







lim i? 2 



T— >+0O 



with 



M = C M e 



-(1+/3)t 



i + o( e - bT ) 



(3.31) 



(3.32) 



where j3 is given by (|3.26|) , Cm is a constant that depends on the initial conditions and b is 
a positive constant. 

Theorems 3.1 - 3.3 enable us to determine the asymptotic behaviour of any physical or 
geometrical quantity in a non-tilted Bianchi VIIo universe, and in particular, provide the 
proofs for Theorems 2.2 - 2.4. The shear parameter S is obtained from equations ( |3.8|) and 
( p.!5| ), and the density parameter Q from equation ( |3.22 ): 



S 2 + R 2 cos 2 ^, n = 1 - s 2 



R 2 



(3.33) 
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Theorems 3.1-3.3, in conjunction with equation ( |3.33 ), now lead directly to Theorems 2.2 
and 2.3. 

Next, the asymptotic form of the Hubble parameter H is determined algebraically, since 
the density parameter Q = fi/(3H 2 ) has a non-zero limit as r — > +oo, for | < 7 < 2, 
and the matter density is given by /x = /io(//^o)~ 37 , or equivalently, /i = /ioe~ 37T , as follows 
from the contracted Bianchi identities. Knowing H, the relation between clock time t and 
dimensionless time r can be obtained by integrating equation ( |1.4| ). On introducing the 
parameter /3, we obtain 



H^H e-^\ ff t « _J_ e CMK (3.34) 

as r — > +00. The asymptotic expansions in Theorems 3.1-3.3 also determine the asymptotic 
form of the line-element by a quadrature and by algebraic means (See Appendix C). 

We now discuss the behaviour of the Weyl curvature. We require an explicit expression 
for the Weyl parameter W, defined by equations fl2.3|) and (|2.4j ). For non-tilted Bianchi VIIo 
models, E ab and H ab , and their dimensionless counterparts 

£ab = j0 , H a b = -jjj, (3.35) 

are diagonal relative to the standard orthonormal frame (see WE, Chapter 6, Appendix). In 
analogy with ( |3.2| ) we define 

= 2 (^ 22 + ^ 33 ) ' £~ = 273 ^ 22 ~~ ^33) 

H + = l(H 2 2 + H 33 ) , H- = ^(H 2 2-H 33 ). (3.36) 

It then follows from equations (|2.3| ) and fl2.4j) that 

W 2 = £ 2 +£ 2 + ?4 + ft 2 . (3.37) 

In terms of the variables R, M,ip as defined by equations ( |3.15[) and (|3.16|) , the 
components £± and TC± are given by 

£+ = E+(l + E+) + |i? 2 (l - 3cos2^) 
ft + = -| J R 2 sin2V' 

£- = jf [sim/> + |M(1 - 2S + )cos^] (3.38) 
H- = ^ [-cos^ - §M£ + sim/>] . 

These expressions follow from equations (6.36) and (6.37) in WE. Since S + and R are 
bounded, it follows from (3.37) and (3.38) that 

W=^[1 + 0(M)], as r^+00. (3.39) 
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This result, in conjunction with Theorems 3.1-3.3, yields 



W 2 



1 



^r- 3 e 2 -[l + 0(^)], if 7 = | (3.40) 



^hfil e 2(l+2/3)r[ 1 + ( e -6r)] j jf | < 7 < 2. 



Theorem 2.4 follows immediately from this result using equation ( |3.27 ). We note that W 
diverges most rapidly for 7 = |. 

As mentioned at the end of Section 2, if | < 7 < 1 the limits of the primary dimensionless 
scalars S 2 , Q and W exist, even though the model is not asymptotically self-similar. We can 
use the asymptotic expressions in Theorem 3.1 to show that if | < 7 < 1, i/iere zs an 
expansion-normalized scalar that is unbounded as r — > +00. We need to know the dominant 
term in the derivatives of and 7i_, where (n) denotes the n th derivative with respect 
to t. One can calculate £_ and 7i_ by differentiating equations (3.38) with respect to r, 
using the evolution equations (3.17) - (3.20). The dominant term arises from differentiating 
the leading trigonometric function, and we obtain 

£ W = f|[cos^ + 0(M)], 

H W = §2 [sin i/) + 0(M)], (3.41) 

as t — > +00. Repeating this process shows that the dominant term in 8^ and grows 
like R/M n+1 as r —>■ +00. Theorem 3.1 implies that 

R 



Mn+l 



„ e Kl-ffl-(2M)] T as r ^ +00 . 



It follows from this result that and are unbounded as r — » +00 if and only if 

n > 6(1 — 7)/(37 — 2). Thus, for any 7 in the range | < 7 < 1, £_ and 7i^^ are unbounded 

as r — >• +00 for n sufficiently large. 

We conclude this Section by giving a heuristic justification of Theorems 3.1 - 3.3. The 
key result underlying the proof is that 

lim M = 0, (3.42) 

as follows from equations ( |3.14j ) and ( |3.16| ). We shall show that in fact M tends to zero at 
an exponential rate. Equation (3.20) implies that once M is close to zero, ip' is large, and 
that ip grows exponentially as r — > +00. Because of this exponential growth, the solutions 
of the evolution equations (3.17) - (3.20) are approximated by solutions of the simplified DE 
obtained by dropping the cos 2ip and sin 2if) terms in equations (3.17) and (3.18): 

t' + = -(2 - Q)t + - R 2 (3.43) 
R! = (Q + S+ - 1)R, 

where we use the symbol * to distinguish the variables in (3.43) from those in equations 
(3.17), (3.18) and ( |3.21| ). This behaviour, which is shown very convincingly in numerical 
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experiments, is due to the fact that the cos 2ip and sin 2ip terms are essentially nullified 
since they change sign repeatedly over an increasingly short time scale. 
The state space for the variables S + and R is 



±1 + R 2 < 1, R > 



see ( |3.24| ) and ( |3.25|) ). Since this set is closed and bounded we can use standard methods to 



analyze the asymptotic behaviour of solutions of the DE (3.43). We find that if | < 7 < ~ 
the point (E+,R) = (0,0) is a global sink for solutions with Q > 0, and if | < 7 < 2, the 
point (£+, R) = \J—(3{1 + (3)) is a global sink for solutions with Q > and R > 0. One 
can also determine the asymptotic form of the solutions of the simplified DE (3.43), which 
gives precisely the leading asymptotic dependence on r in Theorem 3.1 - 3.3. 

In summary, for any initial conditions, equation ( |3.42| ) implies that M will be close to 
zero for r sufficiently large, after which the evolution will be governed asymptotically by the 
simplified DE (3.43). Proving the theorems thus hinges on proving that the oscillatory terms 
do in fact become negligible for r sufficiently large. A detailed proof of Theorem 3.1 is given 
in Appendix B. The proof of Theorem 3.3 ( the case | < 7 < 2) is of a similar nature to that 
of Theorem 3.1, while the proof of Theorem 3.2 (the case of 7 = |) is more complicated, 
since the global sink for the averaged DE (3.43) is non-hyperbolic. The proof thus requires 
the use of cente manifold theory (see for example Carr 1981). The details of this case will 
be published elsewhere. 



4 Discussion 

In this paper we have described for the first time the evolution into the future of any non- 
tilted Bianchi VIIo universe with perfect fluid matter content and equation of state 
p = (7 — l)/i, for 7 in the range | < 7 < 2. We have also given asymptotic forms as r — ► +00 
of various physical and geometrical quantities of interest. Although certain aspects of this 
problem have been studied before, most notably by f\ Collins & Hawking 1973a, Doroshkevich 
et al 1973 and Barrow Sz Sonoda 1986, a general treatment has not been given, and moreover 
a number of important properties of the dynamics have remained unnoticed. 

The most significant feature of the late time evolution is the breaking of asymptotic 
self-similarity (Theorem 2.1), characterized by oscillations that become increasingly rapid 
in terms of dimensionless time r as r — > +00. This oscillatory behaviour leads to the 
phenomenon of Weyl curvature dominance, i.e. expansion-normalized scalars constructed 
from the Weyl tensor become unbounded as r — > +00. We note that this type of asymptotic 
symmetry breaking is quite different from that displayed by the Mixmaster models in the 
singular asymptotic regime, which has a stochastic nature but with all expansion-normalized 
scalars remaining bounded. 

It is worth noting that both types of asymptotic self-similarity breaking create difficulties 
when one attempts to solve the Einstein field equations numerically. In the Bianchi VIIo 
case, if one uses r as the time variable the increasingly rapid oscillations (like cos (e Ar )) cause 

8 We discuss the results of these authors in Appendix C, when we describe the asymptotic form of the 
line- element. 
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numerical difficulties, while if one changes to a time variable relative to which the oscillations 
are asymptotically periodic (e.g. conformal time 77; see table C.I) then the overall evolution 
slows to such an extent that one cannot integrate very far into the late time asymptotic 
regime. On the other hand, as one follows the evolution of a Mixmaster model into the past, 
the model spends an increasingly long time (r-time) in the successive Kasner states so that 
it is likewise difficult to integrate far into the singular asymptotic regime. f\ 

The late time behaviour of the Bianchi VIIo models also has interesting implications as 
regards the question of isotropization. In a spatially homogeneous cosmology, there are two 
physical manifestations of anisotropy, 

i) the shear of the timelike congruence that represents the large scale distribution of the 
matter, 

ii) the Weyl curvature, which can be viewed as describing the intrinsic anisotropy in the 
gravitational field: it determines up to four preferred directions, the principal null 
directions of the gravitational field. 

The shear parameter E and the Weyl parameter W quantify these anisotropies. 

In discussions of isotropization in the literature restrictions are imposed on E, namely 

lim E = 0, (4.1) 

t^+oo 

corresponding to asymptotic isotropization (see for example, Collins and Hawking 1973b page 
324), or 

E < 1, (4.2) 

over some finite time interval 0, corresponding to intermediate isotropization (see Zeldovich 
and Novikov 1983, page 550, and more recently, Wainwright et al 1998, pages 331-2 and 343- 
4). Restrictions are not usually imposed on the Weyl curvature. For the class of non-tilted 
Bianchi models E = characterizes the FL models and hence implies that W = 0, making 
it tempting to conjecture that ( |4.1|) and ([4.2| ) imply the corresponding results for the Weyl 
parameter W. Our analysis shows decisively that this conjecture is false: even though ( |1.1| ) 
holds, the Weyl parameter W can be arbitrarily large for dust (7 = 1) and is unbounded for 
the parameter range 1 < 7 < | (which includes radiation), as r — ► +00. The implication is 
that a model which satisfies ( |4.1| ) or ( |4.2j ) is not necessarily close to isotropy, since W may 
not be small. 

It is of interest to compare the late time evolution of Bianchi VIIo universes with the 
evolution of perturbations of the flat FL model. We can do this by linearizing the Bianchi 
VII evolution equations about the FL model, given by E + = = R. Since the perturbed 
solutions are spatially homogeneous and have zero vorticity, it follows that the scalar and 
vector perturbations are zero, giving purely tensor perturbations (Goode 1989, Section V). 

9 We refer to Berger et al 1997 for a new algorithm for improving the situation somewhat. 
10 For example, the time elapsed since last scattering or the time elapsed since nucleosynthesis. 
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The exact solutions of the linearized equations involve Bessel functions but for our purposes 
it is sufficient to know the asymptotic forms, 



C E e-( 1+ ^ , R lin * Cue-* , M lin * C M e-^\ (4.3) 

Comparing this result with Theorem 3.1, we see that the linearized equations predict the 
correct asymptotic form of R and M, for 7 in the range | < 7 < |. This implies that 
solutions of the linearized equations lead to the same asymptotic form for the Weyl parameter 
W as the exact equations (see equation ( |3.39p ). Thus, the Weyl curvature bifurcation at 7 = 1 
is a property of the linearized solutions as well as the exact solutions. 

On the other hand, it follows from equation ( |4.3|) and Theorem 3.1 that the linearized 
equations predict that S + decays faster than it actually does, i.e. the asymptotic form of 
S + is a non-linear effect. Furthermore, comparing equations (|4.3|) with Theorem 3.2, we see 
that the asymptotic form of all variables in the case of radiation (7 = |) is a non-linear 
phenomenon. 

We now discuss the implications of our results as regards the temperature of the CMBR. 
We have seen that the shear of the Hubble flow, which is described by the two frame com- 
ponents £±, has two noteworthy features, in particular for the case of dust (7 = 1): 

i) the £_ component oscillates increasingly rapidly while tending to zero, as the dimen- 
sionless time r — ► +00, 

ii) the rate of decay of the S + component is a non-linear effect. 
The temperature of the CMBR is given by the formula []] 

r°(l + i: a pK a K^dT , (4.4) 

J Te. 



T Q = T e exp 



where the K a are the direction cosines of a particular null geodesic, = o~ a p/H is the 
dimensionless shear tensor, and r is the dimensionless time. The subscripts o and e refer 
to the present time and the time of emission, respectively. It is thus evident that the link 
between the temperature anistotropy and the shear is complicated, particularly in view of 
the oscillatory nature of the shear, and that determining a reliable bound on the shear 
scalar £ will necessitate integrating the full evolution equations and null geodesic equations 
numerically 0- Without investigating this link further, we can, however, use the fact that 
small £ does not imply small W to make the following assertion: a highly isotropic CMBR 
temperature at the present time does not lead to restrictions on the Weyl parameter. Indeed 
the Bianchi VIIq universes provide an example of cosmological models in which the CMBR 
temperature is highly isotropic but the Weyl parameter is not small. As mentioned in Section 
3, the analysis of Maartens et al 1995b shows that in order to obtain a bound on W one 



11 In papers analysing the anisotropy of the temperature of the CMBR the linearized version of this formula, 
valid for << 1, is usually given. See for example Collins & Hawking 1973a. 

12 The anisotropy in the temperature of the CMBR in Bianchi VIIo models has been investigated by Collins 
and Hawking 1973a and by Doroshkevich et al 1975. In contrast to Collins and Hawking, Doroshkevich et 
al argue that it is not sufficient to use linear perturbations of the FL models. 
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requires an observational bound on the time derivatives of the temperature harmonics, which 
is not attainable in practice. 

We conclude with some remarks interpreting our results from a broader perspective. 
Despite the special nature of the models under consideration - they are spatially homogoneous 
and moreover are not of the most general Bianchi type - the complicated dynamical behaviour 
that we have described may have broad significance, for the following reason. The four- 
dimensional state space of the non-tilted Bianchi VIIo models is an invariant subset of the 
state space for more general models, and by continuity the orbits in an open subset of the 
full state space will "shadow", i.e. approximate, the orbits in the invariant subset. In other 
words, an open subset of more general classes of models will display the features of Bianchi 
VIIo models, at least during a finite time interval. Thus, for example, we can assert that 
there exist Bianchi VIIo models with tilt, and Bianchi VIII or IX models with or without tilt, 
with the property that during some extended epoch the shear parameter S is small but the 
Weyl curvature parameter W is large. More generally, this property will also be displayed by 
spatially inhomogeneous cosmologies, for example G2 cosmologies 0, which contain Bianchi 
VIIo cosmologies as a special case. Based on our experience with the hierarchy of Bianchi 
models, namely that dynamical complexity increases with the dimension of the state space, 
we expect the dynamics of Bianchi models to provide a lower bound, so to speak, for the 
dynamical complexity of spatially inhomogeneous models. Understanding the dynamics of 
Bianchi universes is thus a necessary first step in the study of more general models. 

Appendix A: Proof of Theorem 2.1 

We consider the state space S of the Bianchi VII models defined by the inequalities 

N + > 0, N*- 3N 2 > 0, tt > 0, (A.l) 
(see ( |3l2|) and 0). The function 



Z = \l + ^ +Y Ly with , = |(37-2), (A.2) 



satisfies < Z < +oo on S, and the evolution equations (3.10) - ( |3.11| ) imply that 

Z^_ 4[(S + + t;) 2 + (l-^)S 2 ] 
Z ~ 1 + 

Thus if ~ < 7 < 2, Z is increasing along orbits in S, since < v < 1. 



(A.3) 



13 As an example of the dynamics of a special class of models occurring in a more general class (but 
for the singular asymptotic regime), we note that Weaver et al 1998 have recently shown numerically that 
Mixmaster-like oscillations do in fact occur in G^-cosmologies in the presence of a magnetic field, generalizing 
the behaviour of the magnetic Bianchi VIo cosmologies, whose state space is an invariant subset of the infinite 
dimensional state space of the G2-cosmologies (see Leblanc et al 1995) 
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The set of boundary points of S that are not contained in S is the set S/S, where S is 
the closure of S. It follows from (A.l) that S/S is defined by one or both of the following 
equalities holding: 

n = o , Nl - 3iv! = o. 

By equation (A. 2), Z is defined and equal to zero on S/S. We can now apply the 
Monotonicity Principle (see WE, Theorem 4.12, with "decreasing" replaced by "increasing") 
to conclude that for any point x G S, the w-limit set uj(x) is contained in the subset of S/S 
that satisfies the condition \m\Z(y) ^ 0, where s G S/S and y G S. This subset is the empty 

y-*S 

set, since Z = on S/S, we conclude that u(x) = <fi for all x G S. 
We now prove that 

lim N + = +oo (A.4) 

for each orbit in S. Suppose that (A.4) does not hold. Then there exists a number b > 
such that for any To, there exists a r > tq with N + (r) < b. Since the other variables are 
bounded, this implies that the orbit x(t) has infinitely many points in a compact subset of 
S C M 4 and hence has a limit point in S, contradicting uox) = 0. Thus (A.4) holds. In 
conclusion, we have proved that no orbit in S is future asymptotic to an equilibrium point 
in S and that for every orbit in S, (A.4) holds. 

Appendix B: Proof of Theorem 3.1 

The proof of Theorems 3.1 - 3.3 involves two main steps: 

i) Prove that 

lim (E + ,ii) = { (0 ' 0) ; 

for all initial states. 

ii) Deduce the asymptotic form of (S + , R, M) as r — > +00. 
The evolution equations (3.17) - (3.20) are of the form 

x' = f (x) + fi(x) cos lip 
M' = [m (x) + mxix) cos 2ip + m 2 (x)M sin 2ip]M (B.l) 
$ = ^[l + s(x)M sin V], 

where the arbitrary functions are polynomials in x = (E + ,R). We recall equation ( |3.42| ), 
namely, lim M — 0. 

In the analysis we will encounter scalar functions 

Z = Z{x,M), 



< 7 <| 
< 7 < 2 
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whose evolution equation is of the form 

Z' = [z {x) + zi{x) cos 2ip + z 2 {x)M sin 2vp]Z, (B.2) 

where z , z\ and z 2 are polynomials, and 

Zo (0) = -rj < 0, (B.3) 

It is understood that x, M and ip in (B.2) satisfy the DE (B.l). 
Suppose it is known that all solutions of the DE (B.l) satisfy 

Hi : lim x{r) = 0, 

or the stronger condition 

H 2 : x{t) = 0(e~ bT ), M = 0(e~ br ) as r -> +oo, 

for some constant b > 0. 

These conditions in fact imply that Z tends to zero at an exponential rate. The desired 
conclusions, given in terms of the constant r\ in (B.3), are as follows, 

d : Z = 0{e^ r ' +S)T ) as r -> +oo, for any 5 > 0, 
C 2 : Z = C^e - ^ [1 + 0(e- 6r )] , as r -> +oo, 

where C 2 is a constant that depends on the initial condition. 



Proposition B.l: 



H\ implies C\ and implies C 2 - 



Proof: The idea is to rescale Z so as to " suppress" the rapidly oscillating term z\ (x) cos 2ip 

in (B.2) which may not tend to zero as r — > +00. 

Let 

^= i ^ (B.4) 

1 + \Mzi{x)s\n2^ v 1 

Differentiating (B.4) with respect to r, and using (B.l) and (B.2) leads to an equation of 
the form 

Z' =[z (x)+MB]Z, (B.5) 

where B = B(x,M, cos 2ip, sin 2^) is bounded for sufficiently large r. 

Suppose that H 1 holds. Since M tends to as r — > +00, it follows from (B.3) that given 
5 > 0, 

*b(aO + M5 < -77 + 5, 
for all r sufficiently large, leading to C\. 
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Suppose that H 2 holds. Then (B.5) can be written in the form 

Z' = [- v + 0(e- bT )]Z 

as r — > +00, and C2 follows. □ 

We shall also encounter evolution equations of the form 

Z' = -7]Z + B, (B.6) 

where rj is a constant and B = B(x, M, cos 2ijj, sin 2^) is bounded. By multiplying by 
e -Ti(T-T ) an( ^ integrating from To to r, this DE can be written in the equivalent integral form 

Z( T ) = e-^ T - ro) Z(T ) + I e'^^Bds, (B.7) 

J TO 

where the arguments of B are evaluated at s. In connection with this integral form, we will 
need 

Gronwall's Lemma: 

If v(t) is a non- negative C 1 function, and 

v(t) <C + 5 I v(s)ds, 

J To 

where C, 5 are positive constants, then 

v(t) < Ce s{T - To) for all r > r . 
Proof: See for example, Hirsch & Smale (1974), page 169. 
We can now give the proof of Theorem 3.1. 

Proof of Theorem 3.1: 

We first prove that all solutions of the DE (C.l) satisfy lim x(r) = 0. We write the 

T— + + OO 

evolution equation (|3.11|) for Q in the form 

Q' = 2 [(1 + p)E\ + (3R 2 + R 2 cos 2ip] Q, 



using (|3T21| ) - ( g2|) and ( &2l\ ). 
In analogy with (B.4), we let 



a = — , — % . (b.8) 

l + iM J R 2 sin2^ v ; 



Differentiating with respect to r and using (3.18) - (3.20) leads to an equation of the form 

H' = 2[(1 +f3)T l 2 + + R 2 ((3 + MB)]'U, (B.9) 
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where B is a bounded expression in H + ,R,M, cos2?/> and sin 2ijj, for r sufficiently large. 
Since lim M = and f3 > 0, it follows from (B.9) that there exists tq such that Q > 



T— > + 00 



for all r > ro. Since f2 < 1, equation (B.8) and the bound ( 3.25 ) imply that f2 is bounded 
above. It thus follows that lim f2 exists, say 



r— >+oo 



lim Q = L, (B.10) 

T^+OO 

which, by equations (B.8), and ( |3.42| ) gives lim Q = L. Since Q < 1, we have L < 1. 

Suppose L < 1. Then for r sufficiently large, it follows that + i? 2 = l — Q > |(1 — L), 
Equation (B.9) now implies 

= > 0(1 - L) + 2Si + R 2 MB. 
Since lim M = we have \R 2 MB\ < |(1 — L) for r sufficiently large. It follows that 



n 

n 



= > |/3(1 — L) > 



for r sufficiently large, contradicting (B.10). Thus L — 1, i.e. lim $7 = 1, which implies 



T— >+00 



lim (El + i? 2 ) = 0, i.e. lim x(t) = 0. 

We can now use Proposition B.l, with Hi satisfied, to show that R and M decay expo- 
nentially to zero. The evolution equations (3.18) and (3.19) are of the form (B.2) with the 
constant rj in (B.3) given by 

r] = -2(3 for R 2 , rj = -(1 - (3) for M. 

The proposition thus implies that 

R 2 = ( e (-2 P+ 8)ry M = ( ei - { l-f3 )+S )T^ (Rn) 

as t —>■ +oo. 

The evolution equation (3.17) for S + can be written in the form 

E / + = -(l + /?)fiS + - J R 2 (l + S + )(l- cos 2^) (B.12) 
using (|3.21| ), p.22j ) and ( |3.26| ). Further rearrangement yields 

S' + = -(1 + /3)E+ + B(x, cos 2ip), 

where 

\B\ < (l + P)\Z + \ 3 + Cl R 2 , (B.13) 
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with C\ a constant. By (B.7), this DE has the equivalent integral form 

S + (r) =e- (1+/3)(T - ro) S + (r )+ f e~ il+(3){T - s) Bds (B.14) 

J TO 

where the arguments of B are evaluated at s. By equation (B.ll), 

R 2 < c 2 e^ + ^ for 
Since lim S + = by Hi, we can choose r large enough that 



T— > + 00 



^<7^ for r>r, 
Using these inequalities with (B.13), it follows from equation (B.14) that 
|E + (r)| < Ce(- 2I3+5 ^ + 5 f e (-W+*)<r-) \E+(s)\ds, 
Multiplying by e ( 2 P- 5 ) T yields 

v(t) < C + 5 [ v(s)ds, 



where v(r) = e*- 2/3 ^ r |S + (r)|. Gronwall's lemma thus implies that 

|E+(r)| < Ce { - 2 P +25)r for r > r . 

At this stage we have shown that hypothesis i?2 is satisfied. We can thus use Proposition 
(B.l) to conclude that 

R = C R e- fh [l + 0{e- fiT )], M = C M e- (1 ~ l3)T [l + 0{e- bT )}, (B.15) 

as r — >• +oo (see the discussion leading to equation (B.ll). 

The final step is to deduce the asymptotic form of S + (r). It is necessary to do a change 
of variable in order to suppress the rapidly oscillating term cos 2ip in (B.12). Define 

S+ = S + - |Mi2 2 (l + S+)sin2^. (B.16) 

Differentiating (B.16) with respect to r, and using (B.l) and (B.2) leads to an equation of 
the form 

E' + = -(l + (3)E + -R 2 + B, (B.17) 

where 

B = (l + (3)T,% + (3J: + R 2 + MR 2 B, 



14 



The inequality —(! + /?) < —2/3 + 5 is used to make all constants in the exponentials the same. 
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and B is a complicated expression that is bounded as r — > +00. By (B.7), the DE (B.17) 
has the equivalent integral form 

S + (r) = e^ 1+ ^ T0 ^ + (r ) - f e~^^ R 2 ds - f e-^^Bds. 

J TO J T ( B1 g) 



The rates of decay for £ + , i? and M obtained so far imply that the integral involving R 2 
is the dominant term as r — > +00. Using (B.15), this integral can be evaluated, giving 

[ T e-^ T - s ^R 2 ds = T^e" 2 ^ [1 + 0(e-~ bT )} , (B.19) 

J to H 

as r — > +00, for some positive constant 6. The asymptotic form for E+(r), 

E+ (r) = - I S_ C -^[l + 0(e-*)] > 

as r — > +00, follows from (B.16), (B.18) and (B.19), on noting that the term MR 2 in (B.17) 
is of order e~^ l+ ^ T as r — > +00. The value of the positive constant b in the error term is 
re-defined as necessary. 



Appendix C: The metric approach 

Other investigations of Bianchi VIIo universes have used metric tensor components as basic 
variables. We can obtain the asymptotic behaviour of the metric components directly from 
our results in Section 3, and we now do this for the purpose of comparison. 

For any Bianchi model one can introduce a set of group-invariant and time-independent 
one-forms W a , a = 1, 2, 3, relative to which the line-element has the form 

ds 2 = -dt 2 + g af3 (t)W a W , (C.l) 

where t is clock time along the normal congruence of the group orbits (e.g. WE, page 39). 
For group type VIIo, the one-forms can be chosen to satisfy 

dW 1 = 0, dW 2 = W 3 A W\ dW 3 = W l A W 2 . (C.2) 

If the matter-energy content is a non-tilted perfect fluid, then g a p(t) can be diagonalized, 
and for our purposes it is convenient to use the Misner labelling of the metric components. 

ds z = -dt 2 + l 2 [e-^\W l ) 2 + e 2 ^ ++ ^-\W 2 ) 2 + e 2 ^ + -^-\W 3 ) 2 ] (C.3) 

(Misner 1969, page 1323), where I is the length scale function. The structure relations are 
invariant under a rescaling W 2 — > XW 2 , W 3 — > XW 3 , where A is a constant, which 
implies that there is freedom in redefining I and f3 + according to 

/ -> e 2C l, (3 + -> (3 + + C, (C.4) 

15 This transformation is an element of the automorphism group for Bianchi VIIo (see, for example, Jantzen 
1984). 
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where C is a constant. 

The relations between the metric variables and the orthonormal frame variables are given 
in WE (Chapter 9), for Bianchi models of class A. Specializing equations (10.5) and (10.44) 
in WE to Bianchi VII and using (10.7), (10.11), (6.8) and (6.35) in WE leads to 

N + + v^AL = 1 2(/? + +V3/3~) ; N+ _ = _L e 2(/^-V5/n (C.5) 

HI HI 

Introducing the variables R, M and if> by equations (|3.15|) and (|3.16| ) gives 

tanh(2v / 3/3") = v^^Msin^, (C.6) 



e^ + = ^(l-3/? 2 M 2 sinV). (G.7) 



In addition, WE (10.46) reads 

df3 ± 



->± 



(C. 



dr 

where r is the dimensionless time variable defined by equation ( |1.4j) . We note that the length 
scale I is related to r by equation ( |1.3[) 

We can now use the asymptotic expressions for E + , R, M and if) in Theorems 3.1 - 3.3 to 
find the asymptotic form of (3 + and (3~ , and the relation between clock time t and dimen- 
sionless time r. Since M — > as r — > +oo, equation (C.6) implies that (3~ — > as r — > 0, 
and that for |/?~| < 1 we have /3~ ~ ^RM sinif), which determines (3~ directly. We integrate 
(C.8) to find (3 + up to an additive constant. 

By using the freedom (C.4) in conjunction with equations ([L.3|) , ( |3.34| ) and (C.7) we can 
arrange that the constant Cm in the asymptotic form of M in Section 3 is related to H and 
by 

Cm = H l . (C.9) 
We now give the asymptotic expressions for (3 + and /3~ as r — > +oo in the three cases. 



1 <7<! 



with 



ii) 7 



^ ~ 2(3{l- (3) e ~ WT ' r ~ 5 ( ^ fl oioe- T sin [i/>(r)], (CIO) 



/?+ « -ilnr, (3- « ^H^r^ sin [^(r)], (C.ll) 
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Ill 



with 



V'(t) 



H l T 

In this case the asymptotic form of ip(r) cannot be expressed in terms of elementary 
functions. 

|< 7 <2: 



(3 + « -fir, /T « 7-/9(1 +/9)e- (1+/3)T sin ^(r)], 



(C.12) 



with 



T 



H lo(l + P) 



We now briefly describe the results concerning the late-time asymptotic behaviour of 
Bianchi VIIo universes that have been given in the literature. Comparing the results is made 
difficult by the fact that a variety of different time variables have been used. For convenience 
we summarize the time variables in Table C.l. (see also WE, page 244). The asymptotic 
dependence follows from 



Name 


Definition 


Values 


Dependence on r as r — > +oo 


Clock time t 

BKL T 
Conformal rj 
G 2 -adapted £ 


dt i 

dr H 

dt j3 
dT 1 

— — I 

dr\ 

§ = le- 2 ? + 

d£, 


< t < +oo 
-oo < T < 

< 7] < +00 

< £ < +oo 


£^r/ if | < 7 < |. 



Table C.l Time variables for Bianchi VIIo universes. Dimensionless time r assumes all real 
values, — oo < r < +oo, and t — > +oo in the late time asymptotic regime. The constant (3 is 
given by f3 = ^(4 — 37), where 7 is the equation of state parameter. 



equations ( |3.34| ) and ( |1.3| ) and the fact that (3 + — ► as r — > +00 if | < 7 < | (see (C.10)). 
The £-time was apparently first used for Bianchi cosmologies by Collins & Hawking (1973a, 
see page 328), and subsequently by Siklos (1980). We use the name "GVadapted time" 
because the G3 isometry group admits a two-parameter Abelian subgroup G2, and Ze _2/3+ is 
the length scale in the preferred direction orthogonal to the orbits of the G2. 

Collins and Hawking 1973a analyzed the asymptotic evolution of Bianchi VIIo universes 
assuming that the matter in the model consisted of two non-interacting components, dust 
(7 = 1) and radiation (7 = |). The dust dominates at late times. The tilt was not assumed 
to be zero. They proved that if the initial state is close enough to the Einstein-de Sitter 
model, then the Bianchi VIIo model isotropizes in the sense that a/H — > and (3 — > f3^ as 
t — > +00, where (3^ are constants (see their Theorem 3). Prior to proving this result they 
heuristically determine the asymptotic form of (3 ± , finding that 
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0- « e- a (Acos2£ + £sin2£), + « Ce _Q , 

(see page 332). To facilitate comparison we note that e a = I, H = a, and in their line- 
element in Section IV, A = — 2[3 + , A = 2y / 3/3~, ip = &i = 5 2 = 0. On using Table C.l 
and noting that I = loe T , we see that their result agrees with our result (CIO), when 7 = 1 
i.e. (5 — |. It is of interest that if the field equations are linearized about the Einstein 
de-Sitter model they do not give the correct asymptotic form for (3 + . This fact can be seen 
most clearly from the subsequent paper Collins and Hawking (1973b) in which they give the 
asymptotic form of the solution of the linear field equations for Bianchi VIIo universes with 
dust. In terms of the conformal time r] (see Table C.l), 

/3n = -2/3+ w 4, P22 - /3 33 = 2V3/T w ^-(Bcos2r] + Csm2rj) 
(see pages 316-7). 

Doroshkevich et al 1973 analyzed the non-tilted Bianchi VIIo universes using the line- 
element (C.3) and the BKL time T. For dust (7 = 1), they give the asymptotic form 
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/T « C1T-3 sin (C 2 T-3 + Vo), l 2 e~ w+ « C3T-3 

(see Appendix I, equation (1.20)). To facilitate comparison we note that their metric variables 
7,/i and (AiA 2 )a are related to ours by 7 = Z 6 ,/i = 4^0" and (AiA 2 ) 5 = Z 2 e~ 2/3+ . On using 
Table C.l we see that their result for f3~ agrees with our result (C.10) when 7 = 1 (i.e. 
(3 — |). Their result for [3 + simply confirms that e /3+ — > constant as £ — > +00 but does not 
specify how fast. 

For radiation (7 = |) Doroshkevich et al give the following asymptotic forms 

(3~ » CiTrlsin^(T)], / 2 e"^ + » C 3 ^, (C.13) 



with 



* c 2 ^, 0(t) 



rfT T 2 v 7 /n(f )' 

where Ci,C 2 , C 3 and C 4 are constants (see Appendix I, equations (1.14 - 1.15); we have 
relabelled their constants and have truncated their asymptotic expansions). In order to 
compare with our results we note that 

dip 1 dip 1 - T dip 
~dr = IFWdT ~ I^Hq 6 dT' 

as follows from Table C.l and equation (C.10). In addition Table C.l gives T ~ e~ T , 
suppressing constants, and hence 6 ~ r _1 as r — > +00. It now follows that (C.13) is in 
agreement with our result (C.ll). 

Finally, Barrow and Sonoda 1986, using the equations of Lukash 1975, investigated the 
asymptotic behaviour of non-tilted Bianchi VIIo universes with perfect fluid and 7-law equa- 
tion of state with 1 < 7 < |. They used £-time, up to a factor of 2. By solving the linearized 
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field equations they obtained 

0- ~ ^"T^3(Acos2^ + 5sin2£), 

as £ — > +00, after compensating for the factor of 2 (see their equation (4.144)). This form 
agrees with our result (CIO). In order to facillitate comparison, we note that their \x is 
related to f3~ by /j, = 2\^3(3~, and mention that they redefine \i in their equation (4.138). 
Their parameter a is related to our (3 by a = 2/3/(1 — j3). Barrow and Sonoda argue that 
their linearization procedure is consistent only if 7 < |, suggesting that a bifurcation occurs 
at 7 = |. In contrast, we find that a bifurcation occurs at 7 = |, and that the asymptotic 
form (CIO) is valid for 7 between | and |. 
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